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Abstract 

The formation of beam halos has customarily been described in terms of a particle-core model in 
which the space-charge field of the oscillating core drives particles to large amplitudes. This model 
involves parametric resonance and predicts a hard upper bound to the orbital amplitude of the 
halo particles. We show that the presence of colored noise due to space-charge fluctuations and/or 
machine imperfections can eject particles to much larger amplitudes than would be inferred from 
parametric resonance alone. 
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Beam loss is a major concern for high-current light-ion accelerators such as are needed to 
drive high- intensity spallation neutron sources. Just a tiny impingement, ~1 nA m _1 GeV _1 , 
could generate radioactivation that would preclude routine, hands-on maintenance For 
a 1 mA, 1 GeV light-ion beam, i.e., for baseline beam parameters of the Spallation Neutron 
Source (SNS) presently under construction |2j, this criterion translates to just 1 in 10 6 par- 
ticles lost per meter, a quantity that scales linearly with average beam current. Accordingly, 
a comprehensive understanding of beam-halo formation is imperative. 

Early efforts to identify the fundamental mechanisms of halo formation centered on the 
use of a 'particle-core' model 0,0,^. The basic recognition was that if a uniform- density 
core is made to pulsate, particles that initially lay outside the core and that resonate with 
its pulsations could reach large amplitudes and form a 'halo'. This led to the identification 
of parametric resonance as the essential mechanism of halo formation. A key feature of 
parametric resonance in the context of the particle-core model is a hard upper bound to the 



amplitude that a halo particle can reach Q. Because the particle's orbit al frequency is a 
function of its amplitude, at sufficiently large amplitude the particle falls out of resonance 
with the core and thereby its amplitude ceases from growing further. The prospect that 
the beam halo is 'self-collimating' has led to hope that aperture requirements for beamline 
components might be modest. Smaller apertures are preferred in that, for example, they 
favor higher-efficiency operation of the accelerating cavities. In turn, a wealth of studies 
and a large body of literature has developed over thepast ten years that has centered on 
deciphering the maximum halo amplitude. Refs. Jf| 3,1^,0, ID, U, 12, 3 constitute a small 
sample; Ref. documents a recent halo experiment and alludes to a maximum amplitude, 
though that amplitude could not be measured. 

An urgent question is whether there is any physics not included in the particle-core 
model that could significantly influence the maximum particle amplitude. One feature that 
is unavoidable in real accelerators but is commonly overlooked in simulations is the presence 
of noise. The noise will manifest itself by way of the electromagnetic fields external to 
the beam, which then self-consistently influence the beam's evolving space-charge potential. 
Noise sources could include hardware irregularities that establish fluctuating image-charge 
forces, jitter in power supplies, misalignments and/or asymmetries of beamline components, 
etc. In the context of simulations, it could also include details in the space-charge potential 
that the simulation cannot model precisely. A charged particle will experience all of the noise 
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inherent to the total potential. Moreover, the noise will generally comprise a superposition 
of 'colored' noise, i.e., that for which the autocorrelation time is nonzero. For example, the 
autocorrelation time of noise in the collective space-charge potential could be short, say of 
the order of a plasma period, whereas for hardware irregularities/misalignments it could 
be long, say several betatron (orbital) periods. Herein, by generalizing simple particle-core 
models to include noise, we show that the presence of colored noise can potentially boost 
statistically rare particles to ever-growing amplitudes by continually kicking them back into 
phase with the core oscillation. 

Following the ground-breaking work that introduced the particle-core model , we 
consider particles on radial orbits through an infinitely long, axially symmetric, uniform- 
density beam 'core' that pulsates at a single frequency due to an imbalance, i.e., mismatch, 
between the repulsive, collective space-charge force and the confining external focusing force. 
The core radius oscillates according to an equation of motion for the beam envelope. Upon 
linearizing the envelope equation in terms of the core-oscillation amplitude, one finds the 
solution R(t) — R[l + (M — 1) cos cut], wherein u is the core-oscillation angular frequency 
and M = R(0)/R is the mismatch parameter, i.e., the ratio of initial-to- matched core radii. 
For the uniform- density core, i.e., zero-temperature beam, the core-oscillation frequency is 
to = \^2tt, where f2 denotes the external focusing angular frequency. The particle orbits are 
governed by the dimensionless equation of motion 



the transverse coordinate x is normalized to the radius R of the matched, hence stationary, 
beam; time is multiplied by Q which means all frequencies are expressed as multiples of f2; 
and 6(w) is the Heaviside step function. The second and third terms in square brackets 
govern the motion of the particle when it is inside and outside the core, respectively. This 
model will henceforth be called "Model I" . 

Because Model I is strictly one- dimensional and contains a discontinuity in the form of 
a step function, we shall also study a second model for which the unperturbed beam is a 
spherically symmetric configuration of thermal equilibrium (TE) computed in Ref. Q]. The 
dimensionless equation of motion for this model, henceforth called "Model II" , is 



x + x 1 — 




0; 



(1) 



X = 



f2 2 r 2 + $(r), \l>i = ,u$(ri) smut 
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r = sjx 2 + y 2 + z 2 , n = y/0.8(x 2 + y 2 ) + z 2 ; (2) 

in which the external focusing angular frequency is Q = 1.0001 / \/3. As explained in Ref. Q| , 
the coordinates and time are measured in units of Debye length and inverse plasma angular 
frequency, so the normalization differs from that of Model I. The self-potential $(r) cor- 
responds to "intermediate space charge"; the associated density drops with radius over a 
length scale comparable to that of the quasi-uniform core. The potential \&i is a prolate 
spheroidal perturbation whose strength corresponds to the parameter fi. 

To Models I and II we add fluctuations in the form of Gaussian colored noise such that 
uj — > u(t) = u> + 5u(t), with Su(t) sampling an Ornstein-Uhlenbeck process. Its first two 
moments fully determine the statistical properties of the noise: 

(Su(t)) = 0, (6u(t)6u(ti)) oc exp(-|t - h\/t c ), (3) 

in which t c denotes the autocorrelation time. To keep the models simple, we are choosing to 
add the noise to the core-oscillation frequency; however, we have also confirmed that adding 
colored noise to the external focusing frequency does not significantly change the results. 

After generating a colored- noise signal using an algorithm first presented in Ref. , we 
calculate (|5u;|) which becomes a measure of the noise strength. The influence of noise on halo 
formation should in principle depend on its strength and its autocorrelation time. For two 
choices of autocorrelation time, t c = 1.5r and 12r, r denoting the orbital period of a typical 
halo particle, we investigated a broad range of strengths, specifically 10~ 5 < (\Su\) < 1, with 
the goal of ascertaining to what extent the results may be regarded as generic. Manifestations 
of colored noise that a particle might see are illustrated in Fig. which is provided as an 
aid toward conceptualizing the physical meaning of the noise parameters. Shown there are 
manifestations of noise for a fixed strength (|<5u;|) = 0.01 with t c = 1.5r and 12r, and for 
(\5u\) = 0.1 with t c = 12r. 

In a real beam each individual particle will have its own distinct initial conditions and 
thus experience a manifestation of the noise that differs from that seen by each of the other 
particles. For example, in the axisymmetric Model I, each particle initially occupying a thin 
annulus centered at radius x(0) will experience noise differing from that seen by each of 
the other particles initially in that annulus because the particles start at different angular 
coordinates. The same is true for particles initially occupying a spherical shell centered 
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on radius r(0) in Model II. Accordingly, we adopted a 'survey strategy'. Upon choosing 
initial conditions x(0) and r(0) for Eqs. and (j2J), respectively, and for a specific choice 
of noise parameters, we sequentially computed 10,000 orbits, each experiencing its own 
random manifestation of the colored noise, and we catalogued the maximum amplitudes of 
these orbits. We set the initial conditions of the orbit in Model I at x(0) = 1.20, x(0) = 0, 
and in Model II at r(0) = 1.23, r(0) = 0. In the unperturbed TE sphere of Model II, and 
for realistic proton beam parameters, there are ~ 4 x 10 9 particles per bunch, i.e., ~ 0.6 
nC ^(|. There are ~ 3 x 10 4 particles in the range r = 1.23 ± (0.5 x 10~ 4 ), a thin spherical 
shell centered on r(0) and located well into the Debye tail of the bunch. Accordingly, the 
chosen sample size is realistic. 

For Model I, we examine three values of the mismatch parameter: M = 1.5, 1.3, 1.1. 
Orbits are computed from Eq. first without, then with, the noise of Eq. (jSJ) using a 
variable-time-step integrator. For zero noise, the maximum orbital amplitude |x max | does 
have a hard upper bound in keeping with parametric-resonance arguments, and the upper 
bound depends on the core-oscillation frequency u>o. As Fig. El shows, the particle can reach 
relatively large amplitudes for a wide range of frequencies uq, a consideration that can be 
important in the context of higher-order space-charge modes and harmonics. For nonzero 
colored noise, we present results for which the core-oscillation frequency is fixed at too = \/2, 
the value obtained from the linearized equation for the motion of the beam envelope. By 
design, then, Model I is a direct generalization of the particle-core model introduced in 
Ref. we found that different choices of ujq do not change the essential findings. 

For specified noise parameters, we consider the one particle out of the sample of 10,000 
that reaches the largest amplitude during the integration time of 80r, a time that is repre- 
sentative of the transit time through a 1 GeV proton linac. Results for M = 1.5 and fixed 
t c = 12r are provided in the top panel of Fig. in which |x max | versus (\Suj\) is plotted; 
results for M = 1.3, 1.1 are qualitatively similar. The figure also shows the average \x max \ 
reached by particles in the sample. One sees that over a broad range of noise strengths, 
rare particles are ejected to larger amplitudes relative to the parametric resonance alone. 
For example, a mere 1% fluctuation in the core-oscillation frequency more than doubles the 
maximum amplitude reached compared to the case of zero noise. Interestingly, we found for 
t c = 1.5r that the results are very similar. 

For the analysis of Model II, we fix the perturbation parameter at \i — 0.5. Just as 
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in Model I, with zero noise the particle can reach relatively large amplitudes for a wide 
range of frequencies uio, as Fig. |21 shows. For specified noise parameters, we present in the 
bottom panel of Fig. |H] results for which u = 1.7, a completely arbitrary choice of driving 
frequency; different choices of u) do not change the essential findings. Models I and II 
are distinctly much different, yet the results make clear that the influence of the noise on 
the maximum orbital amplitudes in these models is nearly identical. This is a remarkable 
finding in that we constructed Model II ad hoc, with no predisposition toward matching 
the results of Model I. Accordingly, the influence of colored noise on particle orbits, and in 
particular its role in generating large distended halos in time-dependent potentials, appears 
to be generic. Moreover, the collection of findings suggests that the formation of these halos 
is not particularly sensitive to details in either the governing potential or the noise. 

If the number of particles in the sample is increased with all else being the same, then the 
largest amplitude reached by the single special particle increases. As Fig. |U indicates, once 
the sample size is sufficiently large, the maximum amplitude grows quasi-logarithmically 
with increasing sample size. Only in this restricted sense, and for fixed noise parameters 
and fixed integration time, may it safely be said that there is an upper bound to the halo 
dimension. This particular point has actually been observed, but heretofore unexplained, 
in massive parallelized beam-dynamics simulations of an earlier design of the SNS linear 
accelerator that included a number of machine imperfections j^|. In runs involving 10 4 , 
then 10 5 , then 10 6 , then 10 simulation particles, the maximum extent of the halo increased, 
but it seemed to approach a limiting value with runs above 10 8 particles, a large number 
whose value generally depends on the details of the potential. Inasmuch as these runs were 
self-consistent, the phenomenology they reflect is suggestive of the influence of increasingly 
fine resolution of details in the potential that are beyond the scope of a simple particle-core 
model. They also exemplify that a large number of particles is needed to discern the impact 
of these details on halo formation and structure. 

If the integration time is extended indefinitely, as might be physically representative of 
a storage ring, for example, then there are statistically rare orbits that continue to grow 
to seemingly unlimited amplitudes. Examples of such orbits in Models I and II appear in 
Fig. El These long-time orbits exemplify that there is in principle no upper bound to the 
halo amplitude in the presence of colored noise. 

It remains, of course, to explore further the extent to which this phenomenology applies 
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in real machines. Doing so will involve further simulations of beams in real beamlines; 
as we have seen, machine imperfections will matter. One possibly fruitful approach is to 
extract the coarse-grained, time-dependent potential from the simulations and then add 
noise and pursue a statistical analysis of test particles in parallel to what we have done here. 
Alternatively, the colored noise may be built directly into the simulation itself, although the 
simulation will then need to incorporate a sufficiently large number of particles to garnish 
enough statistics on the halo population. A realistic manifestation of the colored noise would 
need to reflect the machine design, i.e., by properly including imperfections in the fields and 
hardware alignment, and details of the evolving space-charge potential such as a sufficiently 
detailed mode spectrum. Of course, as the beam is accelerated and becomes relativistic, 
space charge and its attendant parametric resonance will become decreasingly important, 
and growth of the halo will thereby be curtailed. 

As a relevant aside, we also analyzed this mechanism in the context of a self-gravitating 
stellar system for which environmental noise from surrounding galaxies will self-consistently 
influence the dynamics. Specifically, we considered a perturbed Plummer model, a config- 
uration for which the unperturbed collective potential scales as (1 + r 2 /3)~ 1//2 |l8j |. and we 
applied the same procedure described herein for Model II. Though it is a restoring force, 
gravity is so weak that, combined with the noise, only a relatively tiny oscillatory perturba- 
tion suffices to pump stars to very large amplitudes. The main point, the generality of which 
is highlighted by the addition of this 'gravitational' example, is that colored noise combined 
with parametric resonance will drive a statistically small number of particles to much larger 
amplitudes than parametric resonance can do on its own. The formation of distended halos 
is thus a general byproduct of collective relaxation of nonequilibrium Coulomb systems. 

This work was supported by the Department of Education under Grant G1A62056. 
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FIG. 1: Example manifestations of colored noise along an orbit for (|<5u;|) = 0.01 and for which 
t c = 1.5r (top) and 12r (center), and for (\Su\) = 0.1 with t c = 12r (bottom). 
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FIG. 2: (top 3 panels) Maximum orbital amplitude vs. core-oscillation frequency loq with zero 
noise in Model I for mismatch parameters M = 1.5, 1.3, 1.1; the horizontal line denotes the initial 
condition x(0) = 1.20 and the vertical line denotes the frequency choice loq = y/2. (bottom panel) 
Same for Model II for perturbation parameter \i = 0.5; the horizontal line denotes the initial 
condition r(0) = 1.23 and the vertical line denotes the frequency choice ujq = 1.7. 
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FIG. 3: (top) For the 10,000-particle sample in the potential of Model I with mismatch parameter 
M = 1.5 and t c = 12r, largest amplitude reached by any particle (triangles) and mean maximum 
amplitude of all particles (diamonds) vs. the dashed line denotes \x max \ for zero noise, 

(bottom) Same for Model II with perturbation parameter \x = 0.5. 
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FIG. 4: Maximum halo dimension vs. sample size N for Model I (top) and Model II (bottom). 
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FIG. 5: Long-time evolution of a large-amplitude orbit given noise with (\5uj\) = 0.01, t c = 12r for 
Model I (top) and Model II (bottom). 
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